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Abstract 



We calculate CP-violating rate asymmetries in the rare radiative decays B — > 
ii'* ± 7 and — > p ±r y. They arise because of the interference between leading-order 
penguin amplitudes and one-gluon corrections with absorptive phases, and provide 
unambiguous evidence for direct CP violation. Complementing earlier studies, we 
also investigate gluon exchange with the 'spectator' quark. The bound state effects 
in the exclusive matrix elements are taken into account by a covariant model, which 
yields a branching ratio BR(B — ► K*j) = (4 — 5) x 10 -5 in good agreement with 
the observed value. The bound state effects increase the CP asymmetry, which is of 
order 1% in the channel B — > K*^ and 15% for B — > fry. 
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1. Introduction 



There is much hope to detect CP violation in B- meson physics |TJ, @, |3fl. It would 
complement existing results in Kaon decays and help solving the origin of this phe- 
nomenon. In fact, the large variety of decay modes of the B- meson makes it possible 
to investigate CP violation from many angles. 

A few decay modes, such as B° — > KsJ/^, allow for a simple, essentially model- 
free prediction of the CP- violating rate asymmetries [[|. These originate mainly 
from the mixing of B° and B°, and are analogous to the e-parameter in K-decay. 
However, one would also like to test direct CP violation in the decay amplitudes, 
such as manifested by the e' parameter (for a recent review of direct CP violation, 
see 0). Theoretical investigations are restricted to inclusive decay modes, which 
are difficult to relate to experimental data or require a calculation of an exclusive 
decay within a phenomenological model. 

It is well known that a non-vanishing asymmetry due to direct CP-violation 
requires two interfering amplitudes with different weak (CKM) and strong (rescat- 
tering) phases. The main difficulty is to calculate the latter. Within the present 
understanding of the basic CP-violating mechanism, consisting in the simultaneous 
contribution of several generations of particles to the physical processes in ques- 
tion, it is natural to locate the strong phases in (penguin) diagrams on the quark 
level, involving virtual heavy and on-shell light particles of all these generations. 
This picture was introduced a long time ago by Bander, Silverman and Soni || 
and underlies most treatments of direct CP violation in charged B-decays; we will 
essentially follow these lines as well. 

A prominent (penguin) decay is B — > K*j, recently observed by the CLEO 
group with a branching ratio of (4.5 ± 1.0 ± 0.9) x 10~ 5 0. With only one hadron 
in the final state it is an interesting candidate: One may hope that the theoretical 
uncertainties of the calculated rate asymmetry 

_ T(B- -> AT*- 7 ) - r(£? + -+ if* + 7) m 
a ° P ~ T(B- -> K*~j) + T(B+ — if*+ 7 ) 1 ' 

are minimal and that the clear signature should make experimental detection feasi- 
ble. 

While the rate asymmetry for the underlying quark process, b — *> 57, vanishes 
up to the one-loop level (as the penguin diagram generates an absorptive phase only 
when the photon is off shell) , Soares || showed that there is a non-zero effect at two 
loops. When true hadronic transitions are considered, also bound state effects must 
be included. They can give rise to additional contributions due to absorptive parts 
of diagrams where additional gluons are emitted from inside of the penguin-loop 
and couple to the other constituents of the hadrons. In this paper we include these 
effects and estimate the rate asymmetry in charged B-decays using a simple model 
for the mesons. 
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2. Perturbative contributions 



We write the effective Hamiltonian for the decay B — > K*^y as 

4G ? 



n. 



-^v c (CM0i c) + C 7 ^)0 7 + ...) (2) 



+ [c-^u). 
where 

Oi q) = (q a YLb Q ) ■ (Sfr/vLqp) (3) 
is the familiar four-Fermi operator and 

°7 = * ^ + m * L ) & ^ ( 4 ) 

is the magnetic dipole operator. L and R are the left and right projectors. The 
operator Oj arises from diagrams with an internal top-quark, and hence its coefficient 
is proportional to v% = VtbV t *. It is convenient to express vt in terms of the two 
independent CKM combinations v,- t = VibV* s (i = u, c) by making use of the unitarity 
relation v t = —v c — v u . This enables us to write the amplitudes for our process in 
the convenient form 

A = v u A u + v c A c . (5) 

For the decay B— > p7 one simply replaces the s-quark in the above expressions for 
the operators and in the CKM factors by the <i-quark. The ellipses in eq. (0) denote 
other operators; in particular, four- Fermi operators whose one-loop matrix elements 
in general contribute to the rate. As pointed out in [[|, their effects can, at least in 
leading logarithmic approximation, be absorbed into an effective coefficient C% (/j). 
For a top-quark mass of 174 GeV, the corresponding values of the Wilson coefficients 
are given by C 2 (5 GeV) = 1.096 and Cy // (5 GeV) = -0.305 [|, H; these val- 



ues include some next-to-leading log (NLL) corrections. Of course, a systematic 
calculation at the NLL level not only requires the knowledge of the coefficients to 
this precision, but also two-loop QCD-corrections to the matrix elements (real and 
imaginary parts) have to be calculated. 

The CP violating asymmetry acp has the form [TT] 



= -Alm[v u vl] \m[A u A* c ] 
aCP \v u A u + v c A c \ 2 + \v*A u + v*A c \ 2 

which implies the well known fact that the absorptive parts of A u and A c must 
be evaluated. These are subleading, i.e. vanish in the leading log approximation 



2 Of course, there are several amplitudes corresponding to the different helicity states of K* and 
7. However, as long as not both the polarization of the photon and the angular distribution of the 
decay products of the K* are observed, there are no other CP-violating observables besides the 
rate asymmetry and we ignore the helicity structure throughout. 

3 The relative sign between C2 and C7 depends on the definition of the covariant derivative; we 
use D = d + ig s T a A a + ieQA 
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and when the matrix element are treated only at tree level. In fact, the imaginary 
parts do not arise from the coefficients Cj, but only from the matrix elements of 
the various operators when evaluated to sufficient order in a s . Therefore, when 
we consider an expansion of the amplitude in a s , it will be consistent to neglect 
subleading terms, except in the imaginary parts in the numerator in eq.(^). In 
particular, these absorptive parts are renormalization scheme independent, which is 
of course not generally true for the real parts of the one-loop matrix elements. The 
leading contribution to the decay rate comes from the tree-level matrix element of 
O7 in Fig. 1 which does not involve the spectator anti-quark. 

At next order in a s there are many terms. In particular, when considering the 
bound state, one would also include, in a perturbative calculation, gluon exchange 
between the 6-quark or its decay products and the spectator line. In fact, in a 



particular treatment of high momentum transfer processes |12| , the gluon exchange 
plays the major role. In this scheme the constituents have essentially no transverse 
momentum which must be provided by the gluon. However, in B-meson decays this 



approach tends to yield too small branching ratios |12j and is likely to serve only 
as corrections to the decay width. To estimate their importance, we have analyzed 
some higher order graphs in the framework of our model (see below) and found that 
they only give corrections of up to order 15% in the rate. We therefore expect that 
an expansion in a s is meaningful and that very soft gluons can be absorbed in a 
suitable wave function. Thus we only include the graph of Fig. 1 to calculate the 
rate. 

However, as pointed out, the higher orders are essential for the asymmetry and 
we must isolate the important ones. At order a° s only the operator O7 contributes 
to the decay (after absorbing the contributions of four-quark operators into C7 ); 
thus, A c and A u are equal and real at lowest order. Including order a s corrections, 
we write A c = A + a s A\ and A u = A + a s v4*. Then, ImfA^*] in eq.(|) reads 

lm[A u A* c ] = a s A° (W* - Im^) . (7) 

This shows explicitly that the internal particles, which generate the absorptive phase, 
must have different masses for a nonzero asymmetry (i.e. the i — u,c quarks from 
must be involved in ImAi). All operators omitted in eq. @ enter with a 
coefficient proportional to v u + v c and therefore do not yield an absorptive phase 
which contributes to the asymmetry at order a s . For the same reason, absorptive 
parts of matrix elements of the operators O7 and Og (the analog to O7, but with a 
gluon instead of the photon) do not contribute; only the loop-level matrix elements 
of O2 with a gluon emitted from the loop are relevant. 

There are basically two classes of diagrams depending on whether the gluon from 
the internal quark loop is attached to the final s-quark (Fig. 2) or to the spectator 
(Fig. 3). 

The first class has already been studied in ref. |J and corresponds to the finite 
part of those QCD radiative corrections which are also responsible for the consid- 



erable operator mixing of O2 into O7 JTOl BJ. For completeness, and ease of use, we 



4 



will recast the result of Soares in the framework of effective operators which also 
gives rather compact expressions. Applying the language of Cutkosky rules, the ab- 
sorptive parts arise from "cuts" which contain the uu or cc pair of O2, respectively. 
Other cut diagrams either vanish or would generate an "universal" phase, i.e. one 
multiplied by v u + v c , which contributes only at higher order to the rate asymmetry. 

The relevant cut diagrams are shown in Fig. 2. They yield local contributions 
proportional to the tree-level matrix element of Oj (i.e. proportional to the contri- 
bution of Fig. 1) 

Im(s-y\Oi l) \b) Fig .2a,b = 9 2 s Qa,bC N p a ,b(mi,m b ) (sj\0 7 \b} tiee , (8) 

where C N = (iV 2 -l)/(2iV) is a color factor (C N = 4/3 for iV = 3 colors). Q a = -1/3 
and Qb = 2/3 are the charges of the quarks to which the photon couples in Fig. 2a,b, 
respectively. The masses of the internal and of the 6-quark are denoted by rrii and 
m b , respectively. To evaluate the factors p a ^ in eq. (§), we apply the Cutkosky rules 
in the form 

(2tt) a r Ni d 3 n , Ni x * 

ImA(6-/) = ^E/n7^t^^(^-E^)A(6-,0A(^/) , (9) 

i n=l \ J n n=l 

where A are the (real) amplitudes of the corresponding subprocesses and the sum- 
mation runs over all possible intermediate states i. 

The diagrams with only the gluon coupled to the quark loop (Fig. 2a) contain 
penguin subdiagrams with insertion of ■ They can be represented by an effective 
sbg vertex 

5= £ t • (10) 

where q is the momentum of the gluon, and V is proportional to the vacuum po- 
larization with internal % = u or c quarks. With this building block the sum of the 
diagrams in Fig. 2a yields 

b t™Ia/ " A 1^2 „2\2 7 2 



where 



p a (mi,m b ) = —— — / ImVl— ) (m b - q ) dq , (11) 



ImV(x) = 2 -^^f-^e(x-4) . (12) 

The sum of the diagrams with the gluon and the photon emitted from the quark 
loop (Fig. 2b) involves an effective sbg'j vertex which can be derived from ref. ||T3| 



1% = -^ry \ie^ a ( q ^Ai 5 +^Ai 6 ) 



+l W qPp * quAi23 + i ^ q%p^Ai 26 . (13) 

p^q 1 p^q 1 J 
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In general, for the gluon and the photon being off shell, the quantities Ai n = 
Ai n (zo, zi, z 2 ) are functions of the three variables Zq = s/m 2 , z\ = q 2 /m 2 and 
Z2 = P 2 /™ 2 , where q and p 7 denote the four- momenta of the gluon and the photon, 
respectively, and s is the invariant mass squared of the internal quark pair. In the 
present situation, where the photon is on shell, the form factors Ai n can be written 
as 

Ai fi (zb,*l,0) = -1 + -^^(Q ( Zo ) - Q Q ( Zl )) ?_(Q_( Zo )_Q_( 2;i )) 

Zq — Zi Zq — Z\ 

z 2 
Ai 6 (z , Zl ,0) = +1 + — (Qo(z ) - Q (zi)) + (Q_(zb)-Q-(zi)) 

Zq — Z\ Zq — Z\ 

Az 23 (^o, zi, 0) = Ai 5 (z o ,zi,0) = -Ai 2 6(z ,zi, 0) . (14) 
The functions Qq and Q- are defined by the integrals 

Q-(x) = / — log (1 — xu(l — u)) , Qo(x) = / du log (1 — xu(l — w)) , 

JO M JO 

(15) 

and have imaginary parts for x > 4. On the sheet Ima; > they are given by 



Img_(o:) = -27rlogf v ^ + ^^ ^ , ] m Q (x) = -J—. (16) 



x 



In terms of these functions, we obtain for the total contribution of the diagrams in 
Fig. 2b 



p b (m u m b ) = — [ dE' s dz- 



EL 



l-E>(l-z)¥ 
x{(2-4^(l-,))^ImQ.(^ 

+ (3^-E>-4E: 2 (l-z))lmg (4)} , (17) 



where E' s mb is the energy of the intermediate (cut) s-quark in the rest frame of the 
decaying 6-quark and z is the cosine of the angle between the intermediate s-quark 
and the photon. In terms of the integration variables the invariant mass squared of 
the internal quark pair reads s = 77t|(1 — 2E' S ). 

At z = — 1 the integral in eq. flT7D is divergent. In an inclusive (free) 6-quark 
decay this divergence would have to be cancelled by including the radiation of addi- 
tional gluons or ss-pairs. In the present exclusive situation this logarithmic diver- 
gence is cut off by requiring that the momentum squared q 2 of the virtual gluon is 
sufficiently off shell in order to avoid double counting with contributions from the 
wave functions. The condition q 2 = —m 2 E' s (l + z) < —A 2 yields the integration 
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intervals 



-1 + 



A 2 






1 


< 


K 


< - 


2m\ 




- 2 


A 2 










< 


z 


< 1 


m\E' s 









We close this section with a physical picture of the diagrams involved. For 
instance, when the internal quark loop contains a cc pair, the above diagrams can 
be viewed as perturbative contributions to the decay of the B-meson into some (sc) 
and (cu) states, e.g. a and a meson. These intermediate hadrons rescatter 
then into the final state K* + 7 by radiating the photon and exchanging gluons 
in all possible ways. This implies the presence of the diagrams (see Fig. 3) of the 
second class, involving the spectator line. They yield non-local contributions and 
their explicit form depends on the model by which the meson states are described; 
the one we use will be presented in the next section. 



3. A model for the exclusive matrix elements 

To describe exclusive decays of B mesons involving a light meson in the final state, 
like B — > -Rev and B —>■ pev, we have developed a covariant model |L4| which 
yields good results for the branching ratio and energy spectrum in these decays; we 
apply it here to the rare decay B — > V'y, where V = K* or p. In fact, since most 
models coincide when the hadronic momentum transfer squared is approximately 
zero, we expect that the results are relatively model independent. We assume that 
both B and K* can be described by two (effective) constituents only: B = (b, q) 
and K* = (s,q). To guarantee covariance of the results we require that the four- 
momenta of the constituents add up to the four-momentum of the bound state. 
Due to binding effects, this condition can only be fulfilled if at least one of the two 
constituents has a variable mass We take the mass m q of the spectator anti- 
quark to be fixed; more specifically we put m q = (for reasons which are related 
to our prescription of the light mesons, see below). Consequently, the 6-quark mass 
becomes momentum dependent 

m b = \jm\- 2p B p q , (19) 

where ps and p q are the four-momenta of the B meson and the spectator, respec- 
tively, and me is the mass of the B meson. Eq. ([19|) only makes sense if the 
momentum \p q \ in the B rest-frame is restricted to \p q \ < itlb/2. As derived in 
detail in ref. [jH|] , the B meson can be represented by a matrix ^> b (in spinor and 
color space) 

^b = C b I ^r^y, L, ^ , M\PB Pq /m B \) S B ® ^=1 N , (20) 
J 2E q {2ny y 2{m B -pBP q ) VN 

4 This is quite similar in the model of Altarelli et al. [|l5| 
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with 

s b = + rn b ) 75 i> q . (21) 
For the wave function 4>b{p) we use a harmonic Ansatz 

<Mp)=exp{V/(2p|)} • (22) 

The normalization factor Cb is chosen such that / d 3 p q \C b& B{\Pq\)\ 2 = (2vr) 3 when 
the integration region is restricted to \p q \ < m B /2. Numerically, Cb ~ \/8 7r 3 / 4 Pp 3 
to a very high precision. The only parameter of the wave function, pp, is fixed in 
such a way that our model yields the right value for the meson decay constant 
Numerical values of pp for different choices of fs are given in table 1. 

While the above construction is adequate for the B meson, it turns out to be 
impossible to describe the light (vector) meson along these lines because whatever 
one chooses for the wave function, it is impossible to get a value for the decay con- 
stant which is compatible with experimental data. (Here we refer to the analogous 
statement given in detail for the pion in ref. fl4j]). Therefore, we use a different 
picture for the light meson, assuming that the meson is represented by two con- 
stituents with parallel momenta p s = ypv and p q — (1 — y)pv, respectively. As the 
vector mesons considered here are much lighter than the B meson, we neglect their 
masses and consequently also those of the constituents. The final state vector meson 
is represented by the matrix 

$ V = C V [ dy<f> v (y)E v (g, -=1 N , (23) 

with 

S y = lim f v ^ v + m v ) , (24) 

where py and ey are the momentum and polarization vector of the meson, respec- 
tively. <pv{y) is the quark distribution amplitude in the vector meson, and the nor- 
malization factor Cy is determined by the value of the decay constant fy, defined 
as (0\u'j fl s\V) = my fye v (fx* ~ f P ~ 216 MeV). Explicitly, we take 

0y(2/) = %(1- 2/) (! + •••) ! C v = fy/(4VN) (25) 

where the ellipses denote deviations from the asymptotic form of the wave function 
and SU (3) Flavour breaking effects (see e.g. ref. [0). 

The exchange of extra gluons may lead to ambiguities (double counting and IR 
divergences) in two situations: In the first, there is a quark line, like the 6-line in 
Fig. 4a, which is off shell because of the exchange of a gluon which can be absorbed 
into one of the meson wave functions. In order to avoid counting this configuration 
a second time in the perturbative treatment, we should impose a condition that the 
quark line is sufficiently off shell. 

In the second, gluons which couple directly to a constituent quark of a meson 
should not be included when they are almost on shell: This case would be degenerate 
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with the contribution from the adequate wave function for a three-particle Fock state 
(e.g. ugs in Fig. 2); since the latter are not included in our simple model treatment, 
we should not include such almost on-shell gluons as well. 

Both conditions will be built into our model by requiring that quarks or gluons 
in the above situations are off shell by at least an amount A. (Of course, the cutoff 
parameter may a priori be different in the two situations; however, for simplicity we 
take the same value A in both cases.) This parameter is part of the definition of 
our model and a value of A = 0.2 — 1 GeV should be physically reasonable. This 
automatically cuts off possible soft and collinear divergencies. 

In the following, we denote by p 7 ipv) an d e 7 (ey) the momentum and polariza- 
tion of the photon (vector meson), respectively. The transition matrix elements of 
the effective operators in eq. @ have the general form (for an on-shell photon) 

(V^f\Oi\B) =em B e» {ie^apPjevPv F i[°i\ + ((Pj e v) Pv^ ~ (p<yPv) ey M ) F 5 [Oi]} 

(26) 

which serves as a definition for the 'form factors' F 15 [Oj]. In the following, when 
the mass of the final-state meson is neglected, we have always i*\ = F 5 due to the 
+ 75) structure of O7; hence, we shall drop the subscript of F (]. 
In the framework of our model for the bound states, the hadronic matix elements 
are determined by 

x < j>J^)fo( y )±Tr\E v M 8b V B Mg 9 ,] , (27) 

y Wig ' iV 

where the trace is in Dirac and color space. The spin projectors and Sy in 
eq. ( p71) are given in eqs. (|21|) and (|2~3j), respectively, and the matrices M s t, and M q / q 
(acting in spinor and color space) are related to the quark-level matrix elements by 

(s(?j\Oi\bq) = u s M sb u b ■ v q M qq ,v q , . (28) 

Note, that in the case of a transition where the spectator is not directly involved, 



M qq i can be written as JT4 



M qq , = 2 lo (2n) 3 5 3 ( Pq -p q ,)® 1 N . (29) 
From the decomposition in eq. (PB|) one readily derives the decay width 



T{B^V 1 ) = ^f m%\F[H eff }\ 2 . (30) 

In the (leading) approximation when only the matrix element of O7 is taken into 
account, one obtains 

T(B^V 1 )^2G 2 F m 5 B a em \v c + v u \ 2 (C e 7 ff ) 2 \F[0 7 ]\ 2 . (31) 



5 For the operator O7, our definition of F coincides with the form factor 2Fi commonly used 
in the literature 



9 



4. Evaluation of F[C>7] and branching ratios 

As mentioned in section 2, only F[0 7 ] contributes to (Vj\TC e ff\B) in a leading-log 
calculation of the decay rate. Evaluating the trace in eq. (|27|) and working out the 
d 3 p q integration by making use of the 3-dimensional delta function (which directly 
relates mj = vn B ^y to the momentum fraction y) we get 

fvm = - Cv J « -JL= Mv) fe(^4^) + o(«.) . (32) 

In order to get a rough idea of the effect of higher order contributions to the 
amplitude for B — > K*j, we consider those order a s corrections to the matix element 
of 07 which involve an extra gluon exchange with the spectator. The corresponding 
diagrams, with the gluon radiated either off the b- or s-quark, are shown in Fig. 4. 
The comparison of the lowest order contributions in our model with these diagrams 
is of particular interest, because the latter would be the leading contributions in the 



approach of ref. ||12|| , where the exchange of hard gluons is essential. 

While the second diagram of Fig. 4 vanishes for m s = 0, the first diagram yields 
the following contribution to -FfC^] 



9 2 s 



F[0 7 ] Fi g A — -^m B C N C B C 
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r , x , N E 2 (m B - 2E g ) (2 - y(l - z)) , 

J N b q 2 J2(m B -E q ) 



where E q is the energy of the spectator in the rest frame f\ of the B-meson, and 
z is the cosine of the angle between the spectator and the A'*-meson. Nt, is the 
denominator of the 6-quark propagator and q 2 is the invariant mass squared of the 
exchanged gluon. In terms of the integration variables we have 

N b = -(l-y)m 2 B + 2m B E q , q 2 = -m B E q (l - y)(l - z) . (34) 



From pure kinematics, both can become zero and generate divergencies in eq. (|33|). 
However, this is avoided by the cutoffs included in our model (see section 3) : When 
the 6-quark propagator, N b , goes to zero the corresponding gluon exchange can be 
viewed as already been taken into account in an adequate wave function of the B 
meson. We may therefore exclude these contributions in the perturbative treatment 
and impose the cutoff |iV b | > A 2 . 

Similarly in q 2 , a vanishing factor (1 — z) (the factor (1 — y)E q is already can- 
celled by the wave function <pv{y) and by E q in the numerator) corresponds to 
the kinematical situation where the gluon is on shell and parallel to the spectator. 
Such kinematical configurations would also enter when we consider higher fock- 
components of the mesons, e.g. ugs. Since we do not go beyond 2-Fock states here, 



Since the model is covariant, the result is of course frame independent. 
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we should also cut off this kinematical region by imposing the (covariant) condition 
\q 2 \ > A 2 . Of course, this leads to a weak (logarithmic) dependence of the results on 
the value of A. We find that for A = 0.2 GeV (1.0 GeV) the non-leading corrections 
of eq. (p3|) contribute at most 23% (6%) to the rate. 

In order to estimate the branching ratio we normalize the rate, evaluated at lead- 
ing order, by the theoretical prediction for the inclusive semileptonic decay width. 
To evaluate the latter, we apply our bound-state model to the decaying B-meson 
and obtain 

G F i t r ,2^2 f d3 Pg i, /I-IM2 m b(Pg) ( m 



V *\ ° B \M\P,\)\ ——TT^ 9[—t-t) , (35) 



31 192tt 3 1 c ° l B J (2tt) 3 ^»\wqu\ mB - Eq ^m b (p q ) 

where g(r) = 1 — 8r 2 + 8r 6 — r 8 — 24r 4 log(r) is the usual phase space function. As 
we work in the leading-log approximation we have neglected QCD corrections in 

The branching ratio is then obtained by multiplying with the measured value for 
BR(B -> Xlv t ) = 0.11: 

T(B -> K*>y) 

BR(B -> Vy) = - x 11% . (36) 

Finally, it is straigthforward to apply our description of the B-meson to the 
inclusive decays B — > X s 7. For the ratio of exclusive to inclusive decay rates, which 
is the more convenient quantity for phenomenological investigations, we obtain 



Rk - = r(B-.*. + 7 ) = <^ 641 |F|Gl1 1 ' (37) 
where m^eff is given by 



m b,eff 



^Cl\M\m 2 ^%rnt{p q ) . (38) 
Neglecting 5 , f/(3)Fiavour breaking effects, R is the same for B — > K* r y and B — ► pj. 



5. Absorptive parts from ImF[02] and rate asymmetries 

For the calculation of the CP asymmetry the imaginary part of FfCV] is needed. 
Corresponding to the contributions from the different classes of diagrams (see Fig- 
ures 2 and 3), we split it into 

F[0 2 ] = F[0 2 ] Ftg . 2a + F[0 2 ] Fig . 2b + F[0 2 ] Fig . 3a + F[0 2 ] Fig . 3b . (39) 

The diagrams in Fig. 2, where no gluon is exchanged with the spectator, are 
readily incorporated into our model. As the effective structure of the relevant sub- 
diagram is proportional to (S7IO7I&) [see eq. (§)], the expression is similar to 
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from eq. 



n 2 



ImF[0 2 )Fig2a,b — "V™ 8 ^ Q a > b ^ N ^ B ^ V 



X 



dy y~ — 4>v{y) 4>b(—^-z — — ) p a ,b{m u m b ) , (40) 



where Q a = — 1/3, Q b = 2/3, and p a ,b have been defined in section 2 [eqs. fllTD and 
([T7|)1. Note that m;, = m By /y depends on y according to our bound-state model; 
and therefore [see eq. ([Of)], Pb is only non-zero for 

A 2 

^<^< X • ( 41 ) 



ni 



B 



In the diagrams with gluon exchange to the spectator, the photon can be emitted 
either from one of the external fermion legs (Fig. 3a) or from the % — u,c quark in 
the fermion loop (Fig. 3b). Diagrams in which the photon is emitted from the b- or 
s-quark line do not develop imaginary parts, because the corresponding momentum 
squared of the gluon, q 2 , is negative. The diagrams in Fig. 3a where the photon is 
emitted from the spectator quark involve the effective sbg vertex of eq. (|T0|) and we 
obtain for their sum 



Q 2 1 f 
lmF[0%>} Fig . 3a = j^-j—Q u C N C B C v dzdE q dy 

J.D7T TRj q J 



E ' , (42) 



2(m B -E q ) 1-2/ V ™l 



where ImV^ is given in eq. (]12[) and the momentum of the gluon is related to the 
integration variables according to 

q 2 = m%(l - y) - 2m B E q + m B E q y(l - z) . (43) 

In the diagrams of Fig. 3b the photon and the gluon is emitted from the fermion 
loop, and the sum of the corresponding two one-particle irreducible subdiagrams is 
represented by the effective sbg'y vertex of eq. (|13|) . After implementing it in the 
bound-state model, we obtain in terms of the form factors Ai n of eq. ([14] ) 

\™F[of] FigM = 9sQuCnCbC v i ^ E 2 ME )M ) 

(>l '" J ./2(m B -E q ) 



[ m B p^q 
+ Ml- 8 )-2 £ ,]( 1+ >, 11 
p^q J q 2 
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where 



q 2 = -m B E q {l - y)(l - z) , 
p y q = ^-[(1 - y)m B - E q (l + z)] . (45) 

The arguments of the functions Ai n (zo, 0) are given by z = (p 7 + q) 2 /m 2 and 
zi = q 2 /m 2 . The cut on the virtuality of the gluon (g 2 < —A 2 ) leads to the restricted 
integration intervals 

A 2 

0< y <!- — 

m% 

A 2 _ m B 

< E„ < 



2m B (l-y) ~ q ~ 2 

A 2 

-1 < z < 1 — r . (46) 

m B E q (l-y) v ; 

We can now evaluate the rate asymmetry acp defined in eq. ([!]). In terms of 
the above form factors F[Ok], the amplitudes A u and A c entering in eq. (^) are 
proportional to 

A t ~C e 7 ff (m B )F[0 7 ]+C 2 (m B )F[0¥] . (47) 
Retaining only the leading (order a s ) terms we get 

ImK^] C 2 (m B ) ImF[Q( M) ]-ImF[Q( c) ] 

acp ~ ~ 2 RT^ • (48) 

To render the dependence on the CKM matrix more transparent it is convenient to 
use the Wolfenstein parametrization |17| which yields 



Im[f u v*] J — rjX 2 (for b s transitions) 

\v u + v c \ 2 ~ \ {1 - P )2 +V 2 (for b -> d transitions) ^ > 

Our numerical calculation of the asymmetries is based on a CKM matrix with 
A = 0.2205 [see eq. and p = -0.3, 0.0, and +0.3 for values of f B = 150 MeV, 
200 MeV and 250 MeV, respectively. We vary rj in the allowed range 0.15 . . .0.5 
(for m t = 174 GeV). For the mass values we use m B = 5.28 GeV and m c = 1.5 
GeV, while m^ff from eq. (^) is given in table 1 for different values of f B . The 
corresponding values of the CP asymmetries, together with the branching ratio and 
R, are given in table 2 for the two choices of the cutoff parameter A = 1 GeV and 
A = 0.2 GeV. 

In our approximation, a s which we choose as a s = 0.2, enters linearly and the 
asymmetry is indirect proportional to (fj) which suffers from a considerable 



uncertainty related to the choice of \x [0, pi. We have also estimated SU (3) Flavour 



breaking effects in the wave functions eq. (|25|) . Using a parametrization as in ref. |L9| . 
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the branching ratio for B — > K* r y (B — > pj) increases (decreases) by up to 30%, 
compared to the values in table 2, while \acp\ decreases (increases) less than 4% 
(7%), respectively. We also note that for the contributions of the diagrams in Fig. 2 
alone, the rate asymmetry is essentially independent of uncertainties from the details 
of the hadronic matrix elements^; this is not the case for the diagrams involving the 
spectator (see Fig. 3). 



6. Discussion 

In this paper we have calculated the rate and the CP-violating asymmetry for 
B— > K*^f. The branching ratio of 4— 5 x 10~ 5 is in good agreement with recent CLEO 



measurements [0] and other theoretical approaches |L8|, [L9|| . With a typical asym- 
metry near 1%, about 10 10 B-mesons are required to resolve acp experimentally. A 
simple rescaling yields for B^ p7 an asymmetry of above 10% and a branching ratio 
of 10 -6 , requiring about 10 9 B-mesons. Similar numbers of B-mesons are needed for 
all such rate asymmetries which test direct CP violation. One expects such samples 



of B-mesons at the planned hadronic B-facilities J2U . 

The asymmetry receives contributions from pure spectator-type decays (consid- 
ered first in ref. ||) and from transitions involving the non-decaying meson con- 
stituents. The two add constructively, a fact stressing the importance of bound 
state effects and possible enhancements of the observables. 

The occurence of both mechanisms is expected if one views the rate asymmetry 
as arising from the rescattering of intermediate DD S or Kir meson pairs. In a quark 
picture, there is no reason to neglect one of the possible gluon exchanges between 
the intermediate c or u quark and the other constituents; taking into account both 
contributions just reflects the symmetry of the problem. If one neglects soft rescat- 
tering of the mesons, the asymmetry is essentially a phase space effect (generated 
by the difference of the u and c quark masses) and probably well represented by the 
quark diagrams. 

The asymmetry depends on a cutoff parameter for the momentum of the ex- 
changed gluon. Choosing it between 200 MeV and 1 GeV leads to a variation of the 
asymmetry by about 20%. We feel that this procedure gives reasonable estimates; 
since a good understanding of the infrared region is still lacking, it is difficult to go 
beyond. 

Within our simple model which takes into account the internal (transverse) mo- 
tion of the quarks, corrections to the branching ratio coming from hard gluon ex- 
change between the quark constituents are small. This confirms the fact that such 
gluons are not sufficient to account for the magnitude of large momentum transfer 
processes (for a discussion of these problems see pi]). 

The rate of B — > K*^ has recently been recalculated, using various QCD sum 



rule techniques [fy], |22], |23| . Our approach yields results most similar to the light- 



7 In our model, the effect of the bound-state wave functions does not completely cancel because 
the argument nib of p in eq. ( ^p|) depends on the quark momenta. 
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^ o Icy 

cone sum rules of ref. |IS]; in particular, the form factor i^O^] scales as m B . In 
contrast, it goes as m B 2 at maximum recoil pi = (tub — my) 2 in ref. |]2~2" |. We 
believe that such a scaling is not correct. Using the Gordon decomposition and the 
known result fs ~ m B 1//2 , one obtains on general grounds that ^[^7] scales as m B 3 ^ 2 
at maximum recoil. 
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Figure captions 

Fig. 1: Leading order contribution from Or, which we use to determine the rate. 

Fig. 2: Order a s matrix elements of O2 with an additional gluon exchange between 
the internal quark loop and the final state s-quark. A cross ( x ) indicates from 
where the photon can be emitted. The dashed line denotes the cut generating 
the absorptive phase. 

Fig. 3: Order a s matrix elements of O2 with an additional gluon exchange between 
the internal quark loop and the spectator quark. A cross (x) indicates from 
where the photon can be emitted. The dashed line denotes the cut generating 
the absorptive phase. 

Fig. 4: Order a s corrections to the matrix element of Or with an additional gluon 
exchange between the b- or s-quark and the spectator. 
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Tables 



f B [MeV] 


p F [MeV] 


m b>e ff [GeV] 


150 


516 


4.68 


200 


643 


4.53 


250 


778 


4.37 



Table 1: Wave function parameter pp and effective 6-quark mass m^// 



fB 

[MeV] 




R 

[%] 


BR 

[io- 5 ] 


a CP for A = 1.0 GeV 

[%] 


a CP for A = 0.2 GeV 

[%] 


B -> K*7 


150 




12 


3.6 


0.23 ...0.76 (0.18 ...0.60) 


0.30 ...1.0 (0.23 ...0.76) 


200 




14 


4.3 


0.24 ...0.81 (0.18 ...0.60) 


0.32 ...1.1 (0.23 ...0.77) 


250 




15 


4.8 


0.25 ...0.85 (0.18 ...0.60) 


0.34 ...1.1 (0.23 ...0.77) 


B -> p7 


150 


-0.3 


12 


0.30 


-2.7 ...-8.1 (-2.2 ...-6.3) 


-3.6 ...-11 (-2.7 ...-8.1) 


200 


0.0 


14 


0.21 


-4.9 ...-13 (-3.6 ...-9.8) 


-6.4 ...-18 (-4.6 ...-13) 


250 


+0.3 


15 


0.12 


-10 ...-24 (-7.2 ...-17 ) 


-14. ...-31 (-9.3 ...-22) 



Table 2: Ratio of exclusive to inclusive decays R, branching ratio BR, and CP- 
asymmetry acp for the decays B — > K*j and B — > pj and for two different cut-off 
parameters A. The numerical values correspond to rj varying in the range 0.15 .. . 0.5 
and a s = 0.2. The weak 77-dependence of the BR(5 — > pj) is not shown. The 
numbers in parentheses correspond to taking into account only the contributions of 
Fig. 2 in F[0 2 }. 
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